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Abstract 

Two points are randomly selected inside a three-dimensional euclidian cube. 
The value I of their separation lies somewhere between zero and the length of a 
diagonal of the cube. The probability density V{1) of the separation is constructed 
analytically. Also a Monte Carlo computer simulation is performed, showing good 
agreement with the formulas obtained. 


1 Introduction 


An important problem in geometry and statistics is: given a convex compact space en¬ 
dowed with a metric, and randomly choosing two points in the space, hnd the probability 
density V{1) that these points have a specified separation 1. The study of this problem has 
a long history |^, and recently gained considerable impetus from researchers in cosmic 
crystallography [p|-|lT9 . 


In a recent paper the functions V{1) corresponding to 2D disks and rectangles were 
The methodology introduced in that work is here extended to a 3D euclidian 


obtained pO 
cube. 


2 Preliminaries 

An euclidian cube with side a is assumed, occupying the location 0 < x,y,z < a in 
a cartesian frame. Randomly choosing two points A and B in the cube, we want the 
probability V{l)dl that the separation between the points lie between I and I + dl. The 
probability density V{1) has to satisfy the normalization condition 

r\/3a 

/ Vil)dl = 1. (1) 

Jo 

The calculation can be shortened if one considers the symmetries of the cube. Really, 
if the points A and B have been chosen, imagine the oriented segment A'B' parallel to 
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AB^ with the tip A' coinciding with the origin O. The other tip B' then lies inside a 
larger cube, with side 2a. Since the probability density V{1) clearly does not depend on 
which octant of the large cube contains B', there is no loose in generality in restricting 
the calculation to the cases where B' is in the octant 0 < x,y, z < a. 

With this assumption, the point B' has cartesian coordinates 

B' = (/ cos 6 cos 0, I cos 6* sin 0, I sin 9 ), 

where both angles 0,0 are bound to the interval [0,7r/2]; here 0 is the azimuthal angle, 
while 6 is the polar angle measured from the z = 0 plane. The corresponding tip B in 
the original segment must lie inside a parallelepiped with sides (see hgure 1) 

/a; := a —/COS0COS0, a —/cos0sin0, lz:=a — lsm9. (2) 


a - I cos6'sin0 



Figure 1 The endpoint B of the segment AB must lie inside the parallelepiped with a corner 
at B'{1, 0, 9). 

The probability V{1, 6, (j))dl dd d(j) that the segment AB has length between I and I + dl, 
azimuth between 0 and 0 + d<f>, and polar angle between 6 and 6 + d6 is then 

V{1, 6, (f))dl dO d(f) = k ly Iz cos 9 dl d9 dcj), (3) 

where A; is a constant and where the assumption O<0, 0<7r/2 stands. Performing the 
angular integrations we shall obtain 

V{1) = JJ V{l,9,(j))d9 d(j), (4) 

and we hnally £x k using the condition (|l|). 

To calculate V{1), three cases need be separately considered, depending on the value 
of I relative to a: namely the cases 0 < / < a, a < / < \/2a, and \/2a < I < \/3a. 

3 The case 0 < / < a 

As is seen in the hgure 2, in this case we effectively have 0mm = 0mm = 0, and (j)max = 
Omax = vr/2. Then 

V{1 < a) = kP [ Iz cos 9 d9 [ Ixlydcj) (5) 
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( 6 ) 


= — 6na^l + 8aP — 

8 

where fc = 8/a® as will be fixed later on. 

z 



Figure 2 The triangular intersection of the 2D sphere with radius I and centre O with the 3D 
cube having side a > 1. 


4 The case a < I < \/2a 

In this case the intersection of the 2D sphere (with radius 1) with the 3D cube (with side 
a) is an hexagonal surface as in figure 3. 


z 



Figure 3 The hexagonal intersection of a 2D sphere with radius I and centre O with a 3D cube 
with a vertex in O and having side a such that a < I < \/2a. 

We note that the arcs of circle drawn on the faces x,y,z = 0 have radius I, while those 
drawn on the faces x,y, z = a have radius \/P — a^. 

For convenience of integration we divide the intersection into two regions. In region I 
we have cos“^(a//) < 6 < sin“^(a/Z) and 0 < </> < n/2. 

In region II we have 9min = 0 and 9max = cos“^(a//). To have (j)min{.9) we note that 
the circle drawn on the face x = a satisfies the equation cos0 cos 6* = a/1, so 

(l)min{9) = cos“^ ( ^ ) =: 4>i{9). (7) 

\ I cos 9 ) 
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On the other hand, the circle drawn on the face y = a satishes sin0 cos 6* = a//, so we 
have 




-1 


I COS 9 


= : (j)2{9). 


( 8 ) 


We then hnd 


V{a < I < \/2a) 


^sin ^{a/l) 

r-wD 

pcos ^{a/l) 


/ cos 9d9 

/ c/0 + 

/ cos 9d9 

/ dcj) 

J cos~^ (a/Z) 

Jo 
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2/^ + — 27ra^(4/ — 3a) — 8a{2f + a^)VP — a? 

+24a^/^ cos“^(a//) , 


where fc = 8/a® as will be hxed later on. 


(9) 


( 10 ) 


5 The case \/2a < I < \/3a 


In this case the 2D sphere with radius I intersects the 3D cube with side a in the triangular 
surface shown in hgure 4. 



Figure 4 The triangular intersection of a 2D sphere with radius I and centre O with a 3D cube 
with a vertex in O and having side a such that \/2o < I < VSa. 


As before, the circles drawn on the faces x,y,z = a have radius The azimuthal 

integration is performed between (j)i{6) and (j)2{9) as in the region II of the preceding 
case, and again sin 6*maa: = o,/l] but now cos 9min = \/2a/l. We then hnd 


^sin l(a/0 f4>2{9) 

V{y2a < I < v3a) = kP 1^ cos 9d9 / lxlvd4> 

^ ’ ycos-l(V2a/Z) hm 

8a{P + a^)\/P — 2a^ — + a^)(/^ + 5a^) + 27ra^(3/^ — 4a/ + 3a^) 

+24a^/sec“^(/^/a^ ~ 1) ~ 24a^(/^ + a^) sec“^ y^/^/a^ — 1 , 


( 11 ) 


( 12 ) 
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where k = 8/a®. 

This value for the constant k derives from the normalization condition (j^), namely, 

[ V{1 < a)dl + [ V{a < I < V2a)dl + [ V{'\/2a < I < '/3a)dl = 1. (13) 

Jo J a J \/2a 

6 Graphs of V{1) 

In figure 5 we present a graph of the dimensionless function aV{l) against the dimension¬ 
less variable I/a. 


a .9(1) 



Figure 5 The probability density V{1) of separation I of pairs of randomly distributed points 
inside a cube with side a. The irregular curve is the output of a corresponding computer 
simulation. 

We note that the function and its first derivative are continuous in the whole interval 0 < 
I < a/So. Nevertheless the second derivative is discontinuous at / = a, as discussed in the 
next section. In the figure a normalized histogram corresponding to 150,000 separations 
between pairs of points randomly selected in the cube is superimposed, for comparison; 
the agreement of the two curves evinces the correctness of the calculation. 


7 Comments 


The integration to hnd V{1) in eqs. (!)-(§) is almost trivial; however, not the same can 
be said about the two other cases, namely in going from to (pT]) and from ([TT|) to (|^. 
A computer assistance appears paramount in these two cases, to conhrm every short step 
in the calculation and simplification of expressions. 

Similarly as in [^, the probability density V{1) and its first derivative are continuous 
throughout the entire range 0 < / < a/So. But the second derivative shows a finite 
discontinuity at I = a, although it is continuous at I = \/2a. 

A remarkable feature of V{1) is its behaviour for large values of /; really, near I = a/So 
we find 


aV{l) = - l/af + 0((y3 - l/af), 

5 


(14) 
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so V{1) is essentially a fifth power of a/3 — I/a. We find that 91% of the separations lie in 
the range I E (0, a), 9% he in the interval I E (a, V^a), and only 0.04% have I > \/2a. 
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